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1. Introduction 



The existence of dynamical symmetries in quantum systems leads to their complete or partial 
integrability and therefore gives the opportunity of a more detailed study of their spectral 
properties. 

We describe here the new method for searching the dynamical symmetry operators based 
on the construction of the isospectral Hamiltonians. The algebraic form of the isospectral 
transformations in Quantum Mechanics is realized by the SUSY algebra M - M. 

The standard one-dimensional Supersymmetrical Quantum Mechanics (SSQM) is gener- 
ated by supercharge operators Q ± , which form the SUSY algebra (together with the Hamil- 
tonian of supersystem H): 

{Q + ,Q-} = H; (1) 

[Q t ,lf\ = 0. (2) 
This algebra is represented by 2 x 2 - matrix supercharges 

Q- = (Q + y=^- o); q ± = Td + W(x) (3) 
and by the superhamiltonian H consisted of two Schrodinger operators: 



Eq.(f|), in terms of components, leads to the intertwining relations for two Hamiltonians: 

fc(V = q + h {2) ; q-h?> = h®q-. (5) 

Therefore the spectra of and almost coincide (up to zero - modes of the operators 
q^) and their wave functions are connected by q ± mappings. 

The extension of one-dimensional SSQM with supercharge operators of higher order in 
derivatives was proposed in [|J, 0. The related SUSY algebra turns out to be a polynomial 
one, namely, the polynomial of the Schrodinger Hamiltonian appears in the right hand side 
ofEq.©. ' ' 

Multidimensional generalizations of SSQM with supercharges linear in derivative are 
known 0,0. In the present paper the two-dimensional SSQM generated by supercharge op- 
erators of second order in derivatives is investigated. In this case, in analogy to 1-dimensional 
SSQM we would expect that the corresponding SUSY algebra becomes polynomial in a su- 
perhamiltonian. However, actually we find that the SUSY algebra (P for a wide set of 
potentials is closed by the diagonal operator of dynamical symmetry (see Sect. 2). It corre- 
sponds to the additional integral of motion for the classical system (Sect. 5). In some cases 
this symmetry operator is of second order in derivatives (Sect. 3) and the system admits the 
i?-separation of variables ||. But, in general it is of 4-th order in derivatives. The particular 
solutions for corresponding potentials and coefficient functions of the symmetry operator 
are displayed in Sect. 4. In Sect. 5 the quasiclassical limit of the superalgebra is studied and 
classical integrals of motion of 4-th order were constructed (up to our knowledge, some of 
them are new (compare to O). 
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2. Two-dimensional Darboux transformations of second 
order in derivatives 



Let us consider the intertwining relations retaining the Planck constant h in the 
Hamiltonian and in the supercharge operators with components of general form: 

q + = ( g -)t = h 2 g ik (x)did k + hCA + B. (6) 

The form of metric gi k (x) is determined by the intertwining relations @: 

dig ik (x) + dig lk {x) + d k g u (x) = 0. (7) 

Its solutions can be easily calculated: 

g n = ax\ + aix 2 + h; (8) 
^22 = 5^1 + 02X1 + 62; (9) 

gu = -\^2ax x x 2 + a x xi + a 2 x 2 ) + b 3 . (10) 

Thus the senior in derivative part of supercharges belongs to the E(2) - universal en- 
veloping algebra ||. We distinguish four different classes in second derivatives: 

qW+ = 1 A + CA + B; (11) 

g (2)+ = a pl + 7 A + CA + B; (12) 

g < 3 >+ = a{J, Pi} + 7 A + CA + B; (13) 

q W+ = a J 2 + (3P 2 + 7 A + CA + B, (14) 

where J and P are generators of rotations and translations, correspondingly, and 

The coefficient functions of supercharge operator (|6|) must satisfy differential equations: 

hdiC k + hd k d + h 2 Ag lk - (V^ - V^)g ik = 0; (15) 
h 2 ACi + 2hdiB + 2hg lk d k V^ - (V^ - V^)Q = 0; (16) 

n 2 AB + tfgikdkdiV® + ncA l v {2) - (yw - y (2) )5 = o. (17) 

From Eqs.(^|) it is evident that the generalized superalgebra gives rise to the symmetry 
operator R for the Hamiltonian H: 

{Q + ,Q-} = R; [R,H] = 0. (18) 

This operator determines a higher-order dynamical symmetry which, in general, can not 
be represented by the polynomial of the Hamiltonian H and of the second-order symmetry 
operator R. As it will be shown in the next Section the latter possibility exists for the 
supercharges ([II]) of the class 1 only. However in all cases the closing of the SUSY algebra 
leads to the integrability of the corresponding dynamical system. 

The general solutions of nonlinear Eqs. flT5|) - ( |I7|) for all four classes of the metric are 
not known. But in some cases the particular solutions can be found in the analytical form. 
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3. Particular solutions of class 1: dynamical symmet- 
ries of second order 



Let us solve the system of differential Eqs. fllBD - ([L7]) for the supercharges with the 
diagonal metric g = diag(—l, —1) (h = 1) : 

q + = -A + CM + B; q~ = (19) 

From Eqs.©, (0) it follows that 

C 2 = (Ci + iC 2 ) 2 = az 2 + 8f3z + r, z = x x + ix 2 ; (20) 

V {2) - B = -a | z | 2 +zfi + zj3 + - | C \ 2 -rj; (21) 

y(2) _ y(i) = dzC + q_Q^ ( 22 ) 

where a, r\ are real constants and (3, 7 may be complex ones. Eq. flTTD can be written in the 
form: 

(C8 Z + C8(g)(S I C | 2 ) = G I C | 2 . (23) 

where 

- OL - - - 

G = a + {d z C){d- z C) --{zC + zC) - 2{(3C + (3C). (24) 
To solve ( p3[) it is useful to change the variables 2 to r!, r 2 : 

f dz r dz r dz r dz , . 

T ' = Jc + Jc- ^ = Jc-Jc- (25) 



Then 



B= j]h' 2 I Glcl2dri+ fci' (26) 



where Ffa) is an arbitrary real function. Thus the irreducible SUSY algebra of second order 
in derivatives is realized for a broad variety of potentials. 

For the class 1 the superalgebra (|i~8l) gives rise to the symmetry operator R of second 
order in derivatives: 

R = H 2 + R + 2 V H, (27) 

where R is diagonal operator: 

( RW \ 

rw ■ 



R 

Its components 



= 2(a I z | 2 +4:((3z + (3z))d z d- z - C 2 d 2 z - C 2 d\ - C{d z C)d z - 

C{d- z C)d- z + {Cd z + Cd z ){B + V {1) ) + B 2 - V {1)2 - 2r]V {1) ; (28) 
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R {2) = 2{a\z\ 2 ^{fiz^r pzyd^-^dl-^dl-Cid^d.- 
Cid-^d-, - (C8 Z + Cd- Z )(B - V {2) ) + B 2 - V^ 2 - 2r]V {2) (29) 

are symmetry operators for the Hamiltonians hP~\ hP\ respectively. 

Let us find explicit expressions for the operators H, Q ± , R. They depend on the values 
of a and (3 in Eq.(20). 



i) If a = 0, (3 0, the constant 7 in (|20"D can be excluded by the translation z 
z — (7/8/5). Then the variables (p5j) read 





and they are related to the conventional M parabolic coordinates. It follows from Eq.( 
that 

2n- \ (3 \ 2 t? -2\ (3 \ 2 t 2 t 2 + F(t 2 ) 



B 



Correspondingly, from Eqs.fl21|), (|22|) we find the potentials V^ 1 ' 2 ': 

yd) _ 1 P ] 2 r t 4 + F(t 2 ) „ 



r 2 + r 2 



(30) 



and from Eqs.(p8|), (p9|) - the components of the symmetry operator R: 



4 I g | 2 
r 2 + r 2 



((-l) l 2<n+ I /3 | 2 t?)t 2 - t(F{t 2 



In terms of variables t%, t 2 the Laplacian 



- (d 2 +d 2 " 



A 



contains the same multiplier as in Eqs.(|30|). Thus the spectral problem for both Hamiltoni- 
ans, h^J = Ety, can be solved by the i?-separation of variables ||. Namely, the corresponding 
eigenfunctions (of and h^) with a given energy E can be expanded into the sum 



* = H Z/ «01n(n)02n(T-2) 



(31) 



where v n are constants, and </>i„(ti), 2n (T 2 ) are solutions of 1-dimensional equations 

I P I 2 [-(£ + ^)r 2 T 2r 1+ | | 2 ^(n) = ^in(ri); 



± 



I /9 I 2 h(J5 + V)r 2 2 + F(r 2 )]0 2n = -^02n(r 2 ). 



(32) 



The upper (lower) sign in Eq.(j32"]) corresponds to h^(h^), and are constants of sepa- 
ration, which play role of the spectral parameters for the symmetry operators Rp-> (R^) 
respectively: 

^0in(Ti)0 2 „(r 2 ) = (A n - 77 2 )0i„(r 1 )0 2n (r 2 ). 



5 



Depending on the properties of the function F(t<i) the degeneracy of energy levels, i.e. 
the dimension of the space {Aj} for a given E can be finite or even infinite. 

ii) Let us describe the second case (5 = 0, a > 0, when the suitable coordinates, 



n = — p= In 

'a 



(33) 



z + Jz* + 1 



r a = — J=ln (34) 



are connected to the conventional elliptic ones. In full analogy with the case i) one can 
obtain that: 

B = 2{fl \ h) ( 2g n/i - \fl - hh + F(t 2 )); (35) 

V {t) = 2(/i * + /g) ((-mdnfl + \fl + ^fa)) - m (36) 



h + h 

/2((-l) < 2a 71 /i + / 1 2 /2)-/ 1 F(7i, 
/i + / 2 



- ^, (37) 



where 



V , r- , , 17 



|2 



/i = - f a exp( v / a • n) H — exp(— y/a • n)J ; (38) 

fz = ^ ^7 exp(zVar 2 ) + 7 exp(-i v / « • r 2 )^ . (39) 

In terms of the coordinates n, r 2 the solution for eigenfunctions can be again decomposed 
into the sum (]3T|) , where now the functions 0i n (ri); 2 „(r 2 ) satisfy the equations: 

+ ~ Ud T1 ft + ^A 2 - (E + rfif^Mn) = ^(n); 

-to) + - (^ + ^)/ 2 ^2n(r 2 ) = -^02n(r 2 ) 

and A^ are eigenvalues of the symmetry operator (j37j). 

The analysis of the case (3 = 0, a < is similar and its result can be formulated as 
follows. It is necessary to replace real (ti,t 2 ) in Eqs.(|33|), (|34"D to imaginary (— ir 2) — iri) 
and, respectively, (/i(ti), / 2 (t 2 ) ) in Eqs. (|38|) , (|39"D to (— / 2 (r 2 ), — /i(tl)). Then the relations 
for 5, V^ (1 ' 2) , # (1 ' 2) are given by Eqs.(g|)-(0). 

Note that for specially chosen a>0,/5 = 7 = the supercharge is factorized into 
the product of two conventional superoperators || that corresponds to the separation of 
variables in the polar coordinates. Thus the supersymmetry, i.e. the intertwinning relations 
between two Hamiltonians, leads inevitably to the hidden dynamical symmetry realized by 
the operator R and moreover to the R - separation of variables in the case when a supercharge 
operator contains senior derivatives in the form of Laplacian. 
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4. Particular solutions of class 2: dynamical symmet- 
ries of fourth order 



I. Let us consider the intertwinning relations (H) for the supercharge (^) with the metric 
g = diag{l, — 1): 

q+ = n 2 {dl-d 2 2 ) + hC k d k + B. (40) 
From Eqs. (|15[)-(P~7D we obtain, that: 

C ± = C 1 tC 2 = C ± (x ± ); (41) 
B = ~(C+C- + F 1 {x + + x-) + F 2 {x + - z_)); (42) 
d_(C_F) = -d + (C + F); F = F 1 (x+ + x_) + F 2 {x + - x_), (43) 

and, respectively, for potentials: 

= ±hc + + c'_) + hcl + c 2 _=) + 

Z o 
-(F 2 (x + - xJ) - Fx(x + + xJ)) + const, (44) 

where x± = xi ± x 2 \ d± = d/dx±. 

Eqs.(|43"D can be solved in certain cases: 

1) Let C_ = 0, then 

C+{x+) = = = ; (45) 

5\ exp(v A ■ x+) + 5 2 exp(— v A ■ x+) 

Fx = o-i<5i exp(v / A • (x- + x+)) + o- 2 <5 2 exp(-v / A ■ (x- + £+)); (46) 

F 2 = a\5 2 exp(\/A • (x + — X-)) + a 2 5x exp(— vA • (x + — (47) 

Here and in what follows the Greek letters stand for constants. Depending on the sign A the 
latter ones may be real or complex. 

2) Let the function F allow the factorization: F = F + (x + ) ■ F_(x_). Then from Eq . ([43|) 
we obtain that: 

x± 

c± = ik ± ~kl F±dx>±i (48) 

and there appear two possibilities: 

a) F ± (x ± ) = e±x±, (49) 

b) F± = a± exp(v / A ■ x±) + 5± exp(— vA ■ x±). (50) 

Note that the constant was omitted in the r.h.s. of Eq. fl49|) since the constant solutions 
F± can be deduced from Eq.(|50"D with A = 0. 
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In the case a) we obtain: 



Fi(x- + x+) = ( ^—(x + + a;_) 2 ; F 2 (x + - xJ) 



6 + e - /_ _ \2. 



and in the case b): 

F x (x + + x-) = o + <7_ exp(v^A ■ {x + + x-)) + 5 + 5- exp(— vA • (x + + a?_)); 
F 2 (x + — X-) = a + 5- exp(\/\ ■ (x + — xJ)) + a_5 + exp(— y/\ ■ (x + — xJ)) 



(51) 
(52) 



In the cases 1), 2a) and 26) potentials can be found by means of Eg . (|44|) . 

For this metric the components of the symmetry operator R (|i~8f ) are given by: 

r^) = iQh 2 dld 2 _ + (T4hd^c_ - c 2 )h 2 d 2 + + (T^d + c + - c 2 + )h 2 d 2 _ + 

2(45 - C + C^)h 2 d + d^ + {Ahd^B - hC+d-C- t BC-)hd + + 

{Ahd + B - hC_d + C + =f BC+)hd_ + 

B 2 + h(C-d + B + C+d-B) + 4h 2 d + d-B. 



(53) 



H. Let us examine now the subclass of degenerate metric g = diag(l,0), when the super- 
charges take the form, 

q+ = h 2 8f + hC k d k + B. (54) 
It follows from Eqs.(|l5D-(|l7D that: 



C x = -x 2 F[ + G x ; C 2 = F 1 
V {1) = h(2G' 1 - t 
K 1 (x 1 )+ K 2 (x 2 ) 



V® = h{2G\ - x 2 F[') + \x 2 2 {F 2 )" - x 2 (F x G x )' + 



V {2) = kx 2 F{' + ~xj{F 2 )" - x 2 {F x G x )' + K 1 {x 1 ) + K 2 (x 2 ); 



B 



~(Gi + x 2 F») + l -G\ - \x\F x F'( + x 2 F x G[ - K x (x x ), 



(55) 

(56) 
(57) 

(58) 



where F X ,G X , K x are arbitrary real functions of variable x x , and K 2 (x 2 ) depends on x 2 only. 

After substitution of Eqs. fl5"5"|) - (|5*5| ) into Eq. (|T7[) we arrive to the equation for the 
coefficient functions: 



h 2 



h 



-^G'l' + -{{G\)" + 2Gf) + G X K[ + 2G' X K X - F X (F X G X )' - G\G\ + 



x 2 



2 



(IV) 



h{F[G'{ + 2G' 1 F 1 ") - G 1 (2G' 1 F 1 / + G X F X 



F[K[ + ~F x (Fl)" - 2G X F X - 2F{'W X + 



+ 



.1 •) 



-G X (F 2 )'" + F[{F X G X )" + 3GiFiFf 



F[(F 2 )"' + F X F" + F x ^(x 2 ) = 0. 



(59) 
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It means that K 2 (x 2 ) is the polynomial of x 2 with constant coefficients: 

K 2 (x 2 ) =m - m x x 2 - \^n 2 x\ - ^m^x\ + ^m±x\. (60) 
Therefore Eq.([59"D is equivalent to four equations: 

-^G- + \{{Glr + 2G^+G l K[ + 

2G\K X - Fi(FiGi)' - G\G\ = miF i; (61) 

2 



fc2 

- h{F[G'l + 2^") - G 1 (2G' 1 F 1 ' + G'/F) - 



F[K[ + ^(F 2 )" - 2G?F X - 2F'{K 1 = m 2 F 1 ; (62) 
^(F 2 )"' + ^(FiGi)" + = m 3 F i; (63) 

iF 1 '(F 1 2 )'" + F 1 F 1 " 2 = m 4 F 1 . (64) 

It happens to be possible to solve the system of differential Eqs. (^1|) - (Q) in particular 
cases. These solutions will be presented for the quasiclassical limit Ti — * in the next Sec. 



5. Quasiclassical limit and integrable systems of class 2 



The quantum dynamical symmetries which were found by the intertwinning method 
in previous Sections have the natural analogs, integrals of motion - in the corresponding 
classical systems. For the class 1 these integrals of motion of 2-nd order are known ||. 

Let us describe certain two-dimensional classical systems of class 2 with additional inte- 
grals of motion of the 4-th order in derivatives. Systems of this type are discussed in the 
book ||, where the list of known particular solutions for the coefficient function of integrals 
of motion was presented. 

I. Let us define momenta p± = —ihd± and take the limit H — ► 0. Then for the Lorentz 
metric we derive from Eqs. 



(]53|) , that the classical Hamiltonian 



+ v c 



cl 1 



where 



Vrl 



\{C\ + Cl) + i(F 2 (x + - x_) - F 1 {x+ + *_)) 



(65) 



reveals the additional integral of motion 



WpIpI + cW + c 2 _pl 



2(F 1 + F 2 )p + p_ + B 2 



(66) 
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The functions C±, Fi, B were obtained for certain cases in Sect. 4. 

H. Let us investigate what kind of integrable systems can be obtained from the SSQM with 
degenerate metric (see Sect. 4). At first, let us define classical functions of supercharge from 
Eq.(HD: 



Qd = ~Pi + iCkPk + B; q d = -p\ - iC k p k + B. (67) 
Thus the Hamiltonian 

h d =Pk + Vd, 

where 

Vd = \xl{Fl)" - x 2 (F 1 G 1 )' + K x { Xl ) + K 2 (x 2 ), (68) 
has the integral of motion 

I = QtiQd = Pi + ( C l ~ 2B )pI + C M + 2CiC 2Pl p 2 + B\ (69) 



The functions C k are defined according to Eq.fl55D and other functions satisfy the system of 

equations (|6l"D - (0) with % — 0. 

Let us display the non-trivial solutions which have been found by now (hi are constants): 
1) Fx — ki ^ 0. Then the constants = = 0, and the function Gx(xx) satisfies 

the equation 

GjdGx , - 

- xi, (70) 



'k - \m 2 G\ 

in turn, the function Kx{x\) can be written in terms of Gx{xx) '■ 

K x {xx) = k 2 G\ 2 + -Gl + k imi Gf J Gxdx'. (71) 

For nonzero values of constants k and m 2 the integral can be represented as a sum of elliptic 
functions of 1-st and 2-nd genus. However in two limits Gx(xx) can be an elementary function 
when a) k > 0, m 2 = 0; b) k = 0, m 2 < 0. 

In the case a), after the redefinition of constants k\, k 2 and an appropriate shift of x 2 , 
we find from Eqs.(B3), fl69|), (fTTD, taking into account (ffOl), that 



T/ klk2 "2/3 , 1 r/ 2 . 3/Cimi 2/3 

Ki = 7—x 2 x x + ~;\k 2 -\ ; \ x x - rrixx 2 + const; (72) 

3 4 k 2 



P z x + klvl + 



2kxk 2 x{ /3 pip 2 + 



kxk 2 2/3 1, , 2 3fcimi 2/ 3 ^1 



(73) 



In the case b) and other solvable cases such as 2) G\ = 0, Fx — k\Xx or Fi = kxx\\ 3) 
Fi = /cia; 2 , Gi = k 2 x\\ 4) iq = k\X\, G\ = k 2 ; 5) Fx = 0, the potentials are obtained as 
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cited in || (i.e. the potentials with separable variables), that is why they are not displayed 
here. 

When comparing the potentials (|65| ) (with the coefficient functions C±, F±, F2 in the cases 
1) and 2b) of Sec. 4) with the list in || we find that they are new and posess the dynamical 
symmetries of 4-th order. 

In conclusion let us formulate the recipe that we have established for the construction of 
integrals of motion in Classical Mechanics. For a given classical Hamiltonian we look 
for the complex function q^(x,p) = (q^)', polynomial in momenta and such that its Poisson 
brackets with hff have the form: 

fe> h< d} = -if{x,P)(Qi to h d ] } = if(x,p)q£, (74) 

where f(x,p) is an arbitrary real function. Ii is evident that then the classical factorizable 
integral of motion I = q^q~i exists. 

This paper was supported by the RFBI grant No. 96-0 1-00535 and by the GRACENAS 
grant No.95-0-6.4-49. 
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